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ABSTRACT 

This  paper  presents  results  that  demonstrate  the  auroral  modeling  capabilities  of  the  Air  Force  Research  Laboratory 
(AFRL)  SAMM2  (SHARC  And  MODTRAN®  Merged  2)  radiance  code.  A  scene  generation  capability  is  obtained  by 
coupling  SAMM2  with  a  recently  developed  Clutter  Region  Atmosphere  and  Scene  Module  (CRASMO),  which 
provides  an  approach  for  rapid  generation  of  time  sequences  and  images  of  radiance  clutter.  Modeled  results  will  be 
compared  to  data  collected  by  the  Midcourse  Space  Experiment  (MSX)1  in  the  IR  and  UV- visible  spectral  regions  during 
an  auroral  event  on  November  10,  1996. 

The  paper  is  organized  as  follows.  We  first  present  a  brief  history  of  the  AFRL  SHARC/SAMM  codes,  leading  up  to  the 
current  version,  SAMM2  v.2.  The  SAMM2  UV-visible  auroral  kinetic  model  will  then  be  described,  followed  by  a 
comparison  of  modeled  results  to  the  MSX  data. 

Keywords:  Aurora,  UV- Visible,  IR,  radiance  scene 

1.  BRIEF  OVERVIEW  OF  THE  SHARC/SAMM  RADIANCE  CODES 

The  current  SAMM2  code  base2  represents  nearly  two  decades  of  active  development  of  an  atmospheric  background 
capability.  Development  began  in  1989  with  the  introduction  of  SHARC,  the  Strategic  High- Altitude  Radiance  code.3 
The  initial  goal  of  this  code  was  to  provide  a  first-principles  model  of  the  high-altitude  non-LTE  (non  Local 
Thermodynamic  Equilibrium)  environment.  SHARC  considered  the  vibrational  states  of  atmospheric  radiators  as  non- 
LTE  and  allowed  user-adjustable  chemical  kinetics  mechanisms  and  rate  constants.  The  original  SHARC  code 
considered  seven  radiators:  three  isotopes  of  C02,  four  isotopes  of  H20,  03,  NO,  OH  and  CH4.  The  code  operated  in  the 
infrared  with  a  spectral  range  of  2-40  pm  with  a  maximum  output  spectral  resolution  of  1  wavenumber.  Interaction  with 
the  code  was  available  through  text-driven  user  interface.  Arbitrary  lines  of  sights  (LOSs)  traversing  altitudes  no  lower 
than  40  km  above  the  ground  were  supported. 

Over  time,  the  capabilities  of  SHARC  were  expanded  to  incorporate  a  number  of  new  phenomenologies.  The  1991 
release  of  SHARC-2  introduced  a  radiance  modeling  capability  for  auroral  conditions.  To  support  this  new  feature, 
important  generators  of  infrared  auroral  radiance  were  introduced,  including  C02(v3),  NO  and  NO+.  Rate  constants  for 
important  reactions  in  auroras  were  obtained  from  the  AFRL  AARC  (Auroral  Atmospheric  Radiance  Code)  program.4  A 
multiple  region  capability  was  implemented  to  allow  for  the  definition  of  the  local,  aurorally  enhanced  region  with  its 
own  atmospheric  chemistry.  SHARC-2  supports  LOSs  traversing  arbitrary  paths  through  the  ambient  and  auroral  regions 
at  altitudes  no  lower  than  40  km  above  the  ground.  The  1993  release  of  SHARC-3  introduced  a  twilight  modeling 
capability,  i.e.,  solar  terminator  condition.  This  was  accomplished  by  introducing  the  concept  of  multiple  atmospheric 
profiles  into  the  code.  Each  profile  can  contain  a  unique  solar  zenith  angle,  allowing  a  smooth  transition  from  day  to 
night.  Along  with  this  addition  was  the  creation  of  an  off-line  tool  capable  of  generating  consistent  atmospheric  profiles: 
the  user-defined  SHARC  Atmospheric  profile  Generator  (SAG)  which  accounts  for  temporal,  geographical  and 
geophysical  variabilities.  The  1997  release  of  SHARC-45  added  an  atmospheric  stochastic  structure  capability  to  model 
radiance  clutter.  With  it,  estimates  of  the  radiance  fluctuations  associated  with  temperature  variations  in  the  atmosphere 
were  obtained. 
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SHARC,  however,  remains  limited  as  a  high-altitude  code  only.  The  SAMM  code  (SHARC  and  MODTRAN®  Merged) 
was  introduced  to  circumvent  this  limitation.  It  integrated  MODTRAN®6,  a  standard  low-altitude  LTE  code,  with 
SHARC  to  produce  a  single  code.  SAMM  maintains  the  flexibility  of  SHARC  while  introducing  the  low-altitude 
capabilities  of  MODTRAN®.  This  new  code  had  a  moderate  spectral  resolution  of  1  cm"1.  Like  SHARC,  it  was  an 
infrared  code  with  a  spectral  range  that  was  extended  down  to  1  micron.  SAMM  is  able  to  support  an  arbitrary  LOS 
from  the  ground  to  the  edge  of  the  atmosphere  will  maintaining  all  of  the  previously  created  capabilities  within  SHARC. 
In  addition,  the  SAG  atmospheric  generator  was  updated  to  provide  atmospheric  properties  for  the  full  altitude  range 
supported  in  SAMM  as  well  as  the  ability  to  generate  MODTRAN®-specific  atmospheres  for  MODTRAN®  calculations. 

The  next  major  revision  to  the  code  base  was  released  in  2003.  SAMM2  was  a  rewrite  of  the  SAMM  code  base  into  a 
single,  seamless  capability  with  a  unified  RT  algorithm.  It  maintained  the  essential  low-altitude  databases  of  the 
MODTRAN®  code  (clouds,  aerosols,  and  surface  reflectance)  as  well  as  the  DISORT  multiple  scattering  capability  and 
MODTRAN®  refractive  geometry  code.  All  of  the  key  features  of  SHARC-4  were  maintained,  including  non-LTE 
chemical  kinetics,  auroral  modeling,  solar  terminator  modeling,  and  stochastic  structure  computations.  SAMM2 
introduced  a  next-generation  full-featured  line-by-line  code  coupled  with  a  novel  treatment  of  line  tails.  As  a  result, 
SAMM2  combines  the  key  capabilities  of  MODTRAN®,  FASCODE,  and  SHARC  and  is  capable  of  producing 
background  scenes  consisting  of  a  large  number  of  unique  LOSs.  At  this  time,  a  new  SAG2  climatology  model  was 
introduced  which  supported  the  NLR  MSISE  2000  model  along  with  additional  water  databases. 

The  next  major  revision  to  the  code  base  is  scheduled  to  be  completed  in  2010.  The  spectral  range  of  all  the 
SHARC/SAMM  codes  up  to  this  point  has  been  the  IR.  For  the  next  revision,  comprehensive  coverage  from  the  infrared 
through  ultraviolet  0.3  to  40  pm  (250  to  33,000  cm'1)  wavelength  region  will  be  provided  for  arbitrary  LOSs  in  the  0-300 
km  altitude  regime.  The  new  code  has  also  optimized  the  usage  of  computational  resources,  resulting  in  a  substantial 
reduction  in  memory  usage  and  a  ten-fold  increase  in  the  computational  speed  of  the  line-by-line  algorithm. 
Additionally,  more  approximate  non-LTE  correlated-^  and  band  models  have  also  been  introduced  to  provide  speedy 
computations  in  applications  where  the  utmost  accuracy  is  not  required. 

It  is  the  extension  of  the  code  spectral  range  into  the  UV- Visible  region  that  will  be  the  focus  of  this  paper.  In  the  next 
section,  we  briefly  describe  the  SAMM2  UV- visible  model. 

2.  SAMM2  UV- VISIBLE  MODEL 

A  UV- Visible  model  was  implemented  in  SAMM2  v.2,  the  most  recent  version  of  the  SAMM2  code,  extending  the 
spectral  range  down  to  0.3  mm.  The  kinetics  for  existing  SAMM2  v.2  IR  species  emitting  in  the  visible  spectral  region 
have  been  augmented  to  support  the  new  visible/ultraviolet  capability.  Existing  SAMM2  v.2  IR  species  (vibrationally 
excited)  emitting  in  the  visible  spectral  region  already  include  LTE  vibrationally  excited  radiators,  such  as  S02  and 
N02,  as  well  as  the  NLTE  vibrationally  excited  radiators,  03,  C02,  H20  and  OH.  Highly  vibrationally  excited  states  of 
water  emitting  in  the  visible  are  formed  by  H20  +  O,  02,  N2.  Additionally,  highly  excited  states  of  OH  (v=6  -  9)  can  be 
formed  by  H  +  03  reactions,  leading  to  the  OH  Meinel  bands. 

A  more  complete  extension  into  UV-visible  region  requires  incorporation  of  NLTE  electronically  excited  states  of  02, 
N2  molecules  and  their  ions,  as  well  as  electronically  excited  states  of  O,  N,  0+,  and  N+.  The  current  version  includes 
electronically  excited  species  and  emission  bands  consistent  with  AURIC  (Atmospheric  Ultraviolet  Radiance  Integrated 
Code)7,  i.e.,  the  N2  electronic  transition  bands  (VK,1PG,  2PG),  N2+  (Meinel  and  1NG)  bands,  02  (ATM,  Herzberg) 
bands,  as  well  as  the  O,  0+  and  N  atomic  lines. 

The  excitation  into  high  electronic  levels  is  triggered  by  solar  radiation  in  the  daytime,  and  by  primary  electron 
bombardment  in  auroras.  The  different  processes  producing  electronically  excited  species,  denoted  by  asterisks,  are 
summarized  as  follows, 

•  Daytime  chemistry  is  driven  by 

-  Photodissociation:  hv  +  02,  N2->  O  +  O*,  N  +  N* 

-  Photoionization:  hv  +  (02  ,N2,0  )->  (02+*,  N2+*,  0+*)  +  e 

-  Photodissociative  ionization: 

•  hv  +  (02,  N2)^(0+*  +  O,  N*  +  N+)  +  e 

-  Electron  Impact  Excitation:  (O,  N2)  +  e  ->  (O*,  N2*)  +  e 

-  Electron  Impact  Dissociation:  N2  +  e  ->N*  +  N**  +  e 


-  Electron  Impact  Ionization: 

•  (02  ,N2,0  )  +e  ^  (02+,  N2+,  0+*)  +  e  +  e’ 

-  Electron  Impact  Dissociative  Ionization: 

•  (02  ,N2)  +  e  ^  (O  +  0+*,  N*  +  N+)  +  e  +  e’ 

•  Under  auroral  conditions,  energy  deposition  by  primary  electrons  ep  replaces  the  solar  flux  hv  in  the  three 
photo-driven  processes  shown  above.  Ionization  and  dissociative  ionization  produce  secondary  electrons,  which 
trigger  ensuing  e-impact  processes. 

•  Resulting  excited  and  ionic  species  from  above  production  trigger  chemical  reactions  leading  to  electronically 
excited  species  that  emit  UV- visible  radiance. 

•  Nighttime  quiescent  radiance,  predominantly  from  oxygen,  is  driven  by  0/02  chemistry: 

-  0+  +  e-»0* 

-  0+  +  02  O  +  02+ 

-  02+  +  e  -»  O  +  O* 

-  o  +  o  +  m^o2*  +  m 

An  example  of  SAMM2  spectral  radiances  in  the  daytime  visible  spectral  range  is  shown  in  Fig.  1. 


Day  time  TH=100km 


Fig.  1.  SAMM2  calculated  daytime  TH=100km  limb  spectral  radiances  of  (a)  OH  Meinel  Band,  02  ATM  band,  (b)  N2+  1NG 
and  Meinel  bands.  01  and  NI  atomic  emissions  are  shown. 


In  the  next  section,  an  application  of  SAMM2  to  model  IR  and  UV-visible  radiances  observed  during  an  aurora  by  the 
MSX  experiments  will  be  described. 


3.  SAMM2  MODELING  OF  MSX  AURORAL  DATA 


3.1  Previous  work 

The  auroral  event  observed  by  the  MSX  satellite  on  November  10,  1996  has  previously  been  modeled  by  several  of  this 
paper’s  authors  Sharma  et  al .8,  as  well  as  by  Strickland  et  al9  Using  a  value  of  -  8  KeV  for  the  primary  electrons 
average  characteristic  energy  determined  by  the  Strickland  group,  we  derived  an  aurorally  dosed  region  that  allowed  us 
to  reproduce  the  Band  A  (6um-llum)  radiances  measured  by  the  SPIRIT  III  Radiometer  instrument  aboard  the  MSX 
satellite.8  The  Band  A  radiance  was  shown  to  arise  from  emission  due  to  vibrationally  and  rotationally  excited  NO 
formed  during  the  aurora.  In  this  paper,  we  have  used  the  same  auroral  region  to  model  the  N2+  1NG  (391.4  nm) 
emission  measured  by  the  MSX  SPIM,  using  the  SAMM2  UV-visible  auroral  model  described  in  section  2.  For 
reference,  we  show  some  of  our  previous  results  from  Sharma.8  The  auroral  region  is  reproduced  below  in  Fig.  2,  and 
the  5.3um  NO  spectra  computed  by  the  SHARC  aurora  model  (now  part  of  SAMM2)  together  with  the  Band  A  filter,  as 
well  as  model  results  and  MSX  data,  are  reproduced  in  Fig.  3. 
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Fig.  2.  Aurorally  dosed  region  [Sharma  et  al.,  (2001)].  Fluxes  are  shown  on  a  latitude/longitude  grid,  characteristic  energy  E0=8KeV. 


Wavelength  (|um) 


Fig.  3.  Aurorally  enhanced  NO  (a)  Emission  from  rotationally  thermal  NO(v=l)  from  O+NO  (b)  Emissions  from 
vibrationally  and  rotationally  hot  NO  from  N(4S,2D)  +  02.(c)  SPIRIT  III  Band  A  data  (d)  SHARC  aurora  model  results 
vs.  SPIRIT  III  data.  From  Sharma,  et  al.8 


In  the  next  subsection,  we  briefly  describe  a  recently  developed  toolkit  that,  coupled  with  SAMM2  results,  allows  for  the 
creation  of  radiance  scenes  that  can  be  compared  to  the  MSX  SPIM  data. 

3.2  SAMM2  scene  generation  capability 


To  generate  radiance  scenes  from  SAMM2  results,  we  have  utilized  a  recently  developed  toolkit  for  SAMM2,  dubbed 
CRASMO  (Clutter  Region  Atmosphere  and  Scene  Module).  CRASMO  defines  the  geophysical  parameters  for  the 
modeled  scene,  generates  the  necessary  profiles  for  the  multiple  atmospheric  regions  including  the  auroral  region, 
incorporates  the  sensor  viewing  angles  to  generate  the  multiple  LOSs  for  SAMM2  runs  and  post-processes  the  outputs  to 
compute  the  scene  radiances.  Examples  of  SAMM2/CRASMO  scenes  are  shown  in  Fig.  4.  Figure  4  illustrates  the  effect 
of  the  primary  electron  energy  on  the  auroral  radiance  limb  profile  in  the  MSX  B2  Band  (-2361  cm"1).  The  larger 
characteristic  energy  caused  a  deeper  penetration  and  a  stronger  emission  from  a  lower  tangent  height. 


(a)  E0=  4KeV  (b)  E0=  8KeV 

Fig.  4.  Effect  of  the  primary  electron  energy  on  the  auroral  radiance  limb  profile  in  the  MSX  B2  Band 
(-2361  cm' 1).  Altitudes  of  Emissions  in  the  B  bands  are  very  sensitive  to  the  Characteristic  Energy  of  the  electrons. 

3.3  Results  from  SAMM2  UV-Visible  aurora  model 

Calculated  results  generated  by  SAMM2  using  input  files  from  the  CRASMO  scene  toolkit  are  shown  in  Fig.  5. 
Comparisons  with  two  selected  SPIM  scenes  show  good  agreement  between  model  and  data. 
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The  occurrence  of  maximum  radiance  at  lower  tangent  height  -  95km  for  time=17:34:24-:29  can  be  explained  by  the 
satellite  moving  closer  to  a  strongly  dosed  region  on  the  near  side  of  the  tangent  point,  i.e.,  the  aurora  moving  further 
away  from  the  tangent  point,  as  shown  for  the  top  two  LOSs  illustrated  in  Fig.  2.  These  results  are  consistent  with 
conclusions  by  Strickland,  et  al.9  and  Sharma,  et  al.8 

SAMM2  spectral  radiances  are  compared  to  the  MSX  SPIM  data9  in  Figs.  6  and  7.  Qualitative  consistency  is  shown 
through  the  UV-visible  spectral  range  from  140  nm  to  580  nm.  The  N2  (VK,  2PG,  LBH)  and  N2+  (1NG)  bands  are  well 
reproduced,  as  well  as  the  O  green  line,  while  the  SAMM2  model  does  not  currently  include  Lyman  (a,P),  OI  (130., 
135,  297.),  0+  and  N+  emissions.  Additionally,  the  auroral  model  does  not  currently  include  N2(X,v)  states  with  v  >  5, 
so  the  LBH  sequence  from  180nm  to  240  nm  is  not  modeled. 


Fig.  6.  Comparison  between  SAMM2  spectral  radiances  and  MSX  SPIM  3  and  4  (see  text). 
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Fig.  7.  Comparison  between  SAMM2  spectral  radiances  and  MSX  SPIM  1  and  2  (see  text). 

4.  CONCLUSIONS 

Ambient  UV- Visible  daytime  and  auroral  chemistry  models  have  recently  been  integrated  into  latest  modern  version  of 
SAMM2,  extending  the  spectral  resolution  to  0.3  um.  Results  from  the  SAMM2  auroral  model  were  previously  validated 
in  the  IR  region  against  the  MSX  NO  data  (Band  A).8  In  this  paper,  SAMM2  aurora  modeling  and  scene  generation 
capabilities  have  been  demonstrated  in  the  UVISI  spectral  region.  Scenes  generated  by  the  SAMM2  auroral  model  of 
N2+  1NG  emissions  are  found  to  be  in  good  agreement  with  the  MSX  391.4  nm  SPIM  data.  Qualitative  consistency  have 
additionally  been  obtained  against  MSX  N2  (VK,  2PG,  LBH)  data. 

The  SAMM2  Aurora  Kinetic  Models  continue  to  be  fine  tuned,  including  incorporating  updated  electron  impact  cross 
sections,  ion  chemistry  reaction  rates,  and  including  a  rotational  NLTE  model  for  NO  radiance. 
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